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ABSTRACT 

We present galaxy cluster mass-richness relations found in the Sloan Digital Sky Survey 
Stripe 82 co-add using clusters found using a Voronoi tessellation cluster finder. These re¬ 
lations were found using stacked weak lensing shear observed in a large sample of galaxy 
clusters. These mass-richness relations are presented for four redshift bins, 0.1 < z < 0.4, 
0.4 < z < 0.7, 0.7 < z < 1.0 and 0.1 < z < 1.0. We describe the sample of galaxy clus¬ 
ters and explain how these clusters were found using a Voronoi tessellation cluster finder. 
We fit an NFW profile to the stacked weak lensing shear signal in redshift and richness 
bins in order to measure virial mass (M 2 oo). We describe several effects that can bias weak 
lensing measurements, including photometric redshift bias, the effect of the central BCG, 
halo miscentering, photometric redshift uncertainty and foreground galaxy contamination. 
We present mass-richness relations using richness measure Nyr with each of these effects 
considered separately as well as considered altogether. We also examine redshift evolution 
of the mass-richness relation. As a result we present measurements of the mass coefficient 
(M 200120 ) and the power law slope ( a) for power law fits to the mass and richness values in 
each of the redshift bins. We find values of the mass coefficient of 8.49 ±0.526,14.1 ± 1.78, 
30.2 ±8.74 and 9.23 ± 0.525 x 10 13 /r*Mg for each of the four redshift bins respectively. 
We find values of the power law slope of 0.905 ±0.0585, 0.948 ±0.100, 1.33 ±0.260 and 
0.883 ±0.0500 respectively. 
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1 INTRODUCTION 

Many properties of a galaxy cluster can be measured directly, in¬ 
cluding cluster richness (the number of galaxies in the cluster), the 
brightness of each of the cluster galaxies and the morphologies 
of cluster galaxies. However some important quantities cannot be 
measured directly but must be inferred from measurable proper¬ 
ties. Mass is one such quantity. Thus in order to constrain mass 
we instead measure quantities that can be related to mass; these 
are called mass proxies. Cluster richness is commonly used as a 
mass proxy, but in order for it to give meaningful results, the re¬ 
lation between mass and richness must be calibrated. This calibra¬ 
tion is referred to as a mass-richness relation. The mass-richness 
relations presented in this paper are given in four redshift bins, 
0.1 < z < 0.4 (low-z), 0.4 < z < 0.7 (mid-z), 0.7 < z < 1.0 (high- 
z) and 0.1 < z < 1.0 (all-z). The main efforts to develop mass¬ 
richness relations previously have focused on clusters at z < 0.4. 

Mass-richness relations are important to galaxy cluster cos- 
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mology as they make it possible to conduct cluster counts as a 
function of mass, allowing c onstraints on cosmological parameters 
such a s Q.m, Os , wo an d de.g. Vikhlinin et alJ2009l : lRozo et al.l 
120101 : IZu et al. 120141) . Higher-redshift (median z ~ 0.6) mass¬ 
richness relations are especially relevant in the context of current 
and future large photometric surveys such as the Dark Energy Sur- 
vey 1 Diehjetaf j20 14) and the Large Synoptic Survey Telescope 
(Abate et al l2012h T These surveys will assemble large samples of 
galaxy clusters at redshifts higher than the average z = 0.25 of pre¬ 
viously found mass-richness relations. Beyond cosmology, mea¬ 
surements of cluster mass at higher redshift will provide knowl¬ 
edge about the evolution of galaxy clusters with redshift and a bout 
the distribution and mass of dark m atter halos de.g. Kauffman et all 
ll995l:lAndreon & Congdodl2014l) . 

Methods of de termining cluste r mass include X-ray 
temperature de.g. Ettori. De Grandi & Molendij |2002h. clus¬ 
ter velocity dispersion ( Allen. Evra rd & Mantz 1201 ll) , 
frequency bias cause d by the Sunyaev-Zeldovich ef¬ 
fect dAde et al.l 120131) and tan gen tial sh ear caused 
by weak lensing dSchneider, Kochanek & Wambgans j 
i2005i) . (For more back ground on lensing theory, see 
especially ISchneider, Kochanek & Wambganssl d2005t ). 
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iNaravan & Bartelmannl d 19971) or iMollerach & Rouletl J2002I) . ) 
Although weak lensing signal can be minimal in single clusters, 
it can be maximized by using stacked shear measurements, 
combining shear signal in bins of similar richness and redshift. 

We observed stacked weak lensing shear as well as clus¬ 
ter richness in the Sloan Digital Sky Survey Stripe 82 co-add 
using a sample o f galaxy clusters found us ing the Voronoi Tes- 
selation method ISoares-Santos et alj 1201 11) . The Stripe 82 co¬ 
add reaches to higher redshift (median sa 0.6) than previous sur¬ 
veys such as the Sloan Digital Sky Survey main sample (median 
z = 0.25). The mass-richness relation with the best statistics to date 
1 Johnston et alj|2007h was based on the maxBCG cluster catalog 
1 Koester et alj2007bD of clusters found in the SDSS main sample. 

In § [2] we summarize previous measurements of mass¬ 
richness relations. In §[3]we describe the properties of our data set. 
In §[4] we describe how stacked weak lensing shear measurements 
were done and how mass was found from shear results. In §[5]we 
describe systematics that introduced uncertainties into our mass 
measurements and how we compensated for them. In this section 
we also describe tests conducted to verify the strength of the tan¬ 
gential shear signal. In §[6]we present our results for mass-richness 
relations for the Stripe 82 co-add including each of the systematics 
individually and then altogether. Finally in §[7]we present an initial 
analysis of the redshift evolution of the mass-richness relation. In 
this paper we take Qm = 0.3, = 0.7, Hq = 100/7 and h = 1. 


2 PREVIOUS MEASUREMENTS 

Ijohnston et al .1 J2007I) developed a mass-richness relation by using 
average shear profiles from stacked weak lensing measurements 
around 130 ,000 galaxy clusters at median redshift of 0.25 found in 
the SDSS djohnston et al.ll2007l ; ISheldon et al.ll2009l) . These clus¬ 
ters were taken from the maxBCG cluster catalog dKoester et alj 
2007b). For richness, they measured A 200 , the n umber of cluster 
galaxies inside radius r 2 oo dHansen et al. 1 12005). To find cluster 
mass they measured stacked weak lensing shear in a set of rich¬ 
ness bins. They then fit a model to these shear profiles. There are 
five terms in their ultimate model of the shear data, including (1) 
the BCG point mass; (2) the NFW profile; (3) the mean surface 
mass profile for miscentered clusters (i.e., those not centered on 
the BCG); (4) the mass of neighboring halos (the two halo term); 
and (5) the non-linear shear effect. After finding cluster masses 
(M 200 ) as a function of richness, Johnston et al. fit a power relation 
to their results, of the form 

M 2 00 =M 200 | 20 (^) (1) 

with a mass coefficient M 2 oo |20 describing the cluster mass at 
A 200 = 20 and a power law slope of a. Their final mass-richness 
relation is 

/w X (1.28±0.04) 

M 2OO = (8.8±1.2xlO 13 /7- 1 M 0 )f^J (2) 

More recently a mass-richness r elatio n was measured in the 
SDSS Stripe 82 co-add dSimet et alj|2012 ). The clusters, being a 
subset of the clusters used in Johnston et al., have the same me¬ 
dian redshift of z ~ 0.25, but many more source galaxies can be 
observed as the Stripe 82 co-add can reach magnitude 23 in i-band 
(with 50% completeness). Simet et al. consider several system¬ 
atics, including error introduced by treating the brightest cluster 
galaxy (BCG) as being in the center of the halo when it is not. 


This is called halo miscentering. They consider miscentering of 
the BCG by creating a set of mock catalogs from their data. These 
catalogs preserve galaxy positions, shape errors, photometric red- 
shifts and more but replace actual source galaxy shears with ex¬ 
pected shears from a shear model. Since the shear model is pro¬ 
duced by a halo model, they control how many of the halos have 
a miscentered BCG and can obtain information on the actual halo 
masses. They then fit an NFW profile to the mock data and find 
that the measured mass is underestimated due to miscentering of 
the BCG. They then find a relationship between the richness and 
the scale of miscentering: 


MlOO,true 
T/ooO.m/v 


1.44T0.17 



-0.21±0.18 


(3) 


They also find a mass-richness relation of 


M 200 = (9.56 ±0.75 x 1O 13 M 0 ) 


^y iio±oi2) 


(4) 


iFord et all l l2014al ) describe measurements of lensing mag¬ 
nification in the Canada-France-Hawaii Telescope Lensing Sur¬ 
vey (CFHTLenS) using cluster candid ates fo und using the 3D- 
Matched-Filter cluster-finder of MilkeraidxeUth d2010l) . In an¬ 
other analysis of the same data lFord et alJ J20141: ) describe mea¬ 
surements of weak lensing shear in the CFHTLenS. Measurements 
are made for more than 18,000 cluster candidates at 0.2 < z < 0.9. 
In the first paper the authors fit NFW profiles to the lensing magni¬ 
fication measurements and use these to find best-fitting values for 
M 2 oo. These values are then used to produce a mass-richness re¬ 
lation. Systematics considered include cluster miscentering, pho¬ 
tometric redshift errors and the two-halo term (that is dark matter 
structures nearby the cluster). The final mass-richness relation ob¬ 
tained is 


M 20 o = (2.2 ±0.2 x 10 13 Mq) 


A 2 qo \ 

20 J 


(1.5±0.1) 


(5) 


In the second paper, the authors stack weak lensing shear signal, 
fit the result to an NFW profile and correct for miscentering. They 
find a final mass-richness relation of 


M 200 = (2.7l°; 3 x 1O 1J M 0 ) 


vl3, 


A 2 qo \ 

20 J 


(1.4±0.1) 


( 6 ) 


iReves et aD J2008I) describe a study of ~ 13,000 clusters 
from the maxBCG catalog at 0.1 < z < 0 .3. While this is the same 
sample as that presented bv ljohnston et al J l i2007 :). their approach 
is different in that they try different combinations of mass trac¬ 
ers to identify the optimal mass tracer for galaxy clusters. They 
find cluster masses using stacked weak lensing shear measure¬ 
ments. They consider systematics including photometric redshift 
errors, cluster miscentering, contamination of shear signal and in¬ 
trinsic alignment of source galaxies. They consider the reliability 
of several mass proxies, including richness (Niqo), cluster lumi¬ 
nosity (±200 ) and BCG luminosity (±gcc)’ ultimately combining 
N 20 Q and Lbcg into a mass-proxy relation. Along the way they 
provide a mass-richness relation of 

r N \ (1.164=0.09) 

M 200 = (L42± 0.08 xlO 14 M 0 )(^J (7) 

In Section [6.2.21 we compare our results to the previous results 
discussed in this section. 
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3 DATA 
3.1 Stripe 82 

The Sloan Digital Sky Survey (SDSS) began in 1998 and sought 
to image 10,000 deg 2 in the North Galactic Cap and to take spec¬ 
troscopy of one million galaxies a nd one hundred thousand quasars 
in this same region lAnnis etalj|2014h . The SDSS uses a 2.5-m 
telescope located at Apache Point Observatory. The SDSS camera 
has 24 2048x2048 CCDs of pixel scale 0.396 arcsec. Stripe 82 is 
a section of the footprint of the Sloan Digital Sky Survey (SDSS); 
Stripe 82 falls along the celestial equator between declinations of 
— 1.25° and +1.25°. Stripe 82 was imaged repeatedly during the 
fall when the North Galactic Cap was not observable. The main 
SDSS area reaches r ~ 22.4 and median seeing of 1.4 arcsec, while 
Stripe 82 reaches about 2 magnitudes deeper and achieved median 
seeing of 1.1 arcsec. This is possible as about 20 runs were taken 
of each field in the Stripe 82 co-add and these data were then co¬ 
added. 


3.2 Cluster Samples 

The maxBCG cluster catalog dKoester et al.ll2007al) is a catalog of 
13 823 clusters found in the Sloan Di gital Sky Survey main survey 
using the maxBCG method dKoester et alj|2007bi) . The maxBCG 
method involves searching for clusters based on the presence of 
a brightest cluster galaxy (BCG) and a number of E/S0 ridgeline 
galaxies nearby of similar color. The maxBCG clusters were found 
in a region of area 7500 deg 2 and have 0.1 < z < 0.3. Values of 
N 200 range from 10-190. We used the maxBCG catalog as a control 
sample at low z, verifying that our measurements of weak lensing 
shear reproduced the mass-jdchness relati on pre viously found from 
the maxBCG full sample (Johnston et aljl2007t) . We also used the 
maxBCG cluster sample to find a relation between /V 200 and our 
richness measure Nyj and we utilized a correction for halo mis- 
centering based on maxBCG cluster data. 

We produced a catalog of 19 316 clusters from the Stripe 
82 co-add using a Vor onoi tessellation (VT) cluster finder 
dSoares-Santos et al.l201 ll) . A sample of four higher-richness clus¬ 
ters from the sample are shown in Figure Q] A description of the 
devel opment and the properties of this cluster catalog is in prepara¬ 
tion dSoares-Santos et aill20 1 5l) . Voronoi tessellation is a method 
of relating distances between different locations in a large array 
of points (called seeds). In a Voronoi tessellation, regions called 
Voronoi cells are defined in which all points within that cell are 
closer to that seed than to any other seed. In this method, the seeds 
would be galaxies. If the number and size of Voronoi cells exceed 
certain criteria, then this region is declared to be a cluster. This was 
found to be a way of assembling a cluster sample with high pu¬ 
rity and completeness while not de pending on galaxy magnitude 
or color dSoares-Santos et ail 1 201 ll) . The total number of cluster 
members is defined by the number of objects within the overdense 
region in the VT method. Thus the richness measure is not N 200 
but is instead a new richness measure called Nyj. 

We wanted to find an approximate relation between /V 200 and 
Nvt both to compare results for each and in order to use Equation 
[3] for miscentering corrections (this equation depends on V 200 )• 
In order to do this, we found galaxy clusters listed both in the 
maxBCG catalog and in the Stripe 82 VT catalog. We did this by 
searching the maxBCG catalog for clusters within 1.0 arcmin of 
the cluster center in the VT catalog. We found 230 galaxy clusters 
in both catalogs using this method. We plotted the N 200 from the 



Figure 1. Images of the central regions of four of the higher-richness 
clusters from the new Voronoi tessellation cluster sample in the Stripe 
82 co-add. These clusters (proceeding clockwise from the top left) 
are at (319.70833. 0.55727617), (328.61234, 0.64274948). (16.227354, 
0.066838292), and (345.63669, 0.040305327). They have N V t of 99, 
78, 127 and 76 and z of 0.28, 0.33, 0.30 and 0.69 respectively. Im¬ 
ages are taken from the SDSS DR7 Stripe 82 Data Archive Server at 
http://das.sdss.org/www/html/imaging/dr-T5.html Each im¬ 
age has size of approximately 2.5 arcmin square. 


Comparison of N w and N 200 



Figure 2. A plot of Nvt vs. N 200 for 230 galaxy clusters matched between 
the maxBCG cluster sample and our VT cluster sample. 


maxBCG catalog against the Nyj from the VT catalog; the result 
is shown in Figure [2] There is not a clear relation between Nyj 
and A^oo- Nevertheless we tried to find an approximate relation by 
fitting a power law to the data (the power relation was a better fit 
than a simple linear fit). The equation we obtained relating the two 
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cuts were made on magnitude, size and ellipticity components 
llin et alj2012h . There are 5875 133 objects in this catalog, rang¬ 
ing from 0.3 < z p hot < 1-60 and i-band magnitude 18 < i < 24. 
All galaxies in this catalog include shape information measured 
by the SDSS pipeline. This catalog also contains photometric red- 
shifts measured for each galaxy using the artificial neural network 
technique dReis et al J|20 I j) . These photo-zs were trained using a 
sample of about 83 000 galaxies that have spectroscopic redshifts 
measured for them. 


4 SHEAR MEASUREMENTS 


Figure 3. Histograms of richness for the maxBCG and the Stripe 82 VT 
cluster catalogs. Histograms are normalized by their respective maximum 
values. 
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Figure 4. Histograms of redshift for the maxBCG and the Stripe 82 VT 
cluster catalogs. Histograms are normalized by their respective maximum 
values. The peak at high redshift is due to limitations in the photometric 
redshifts of the input galaxy catalog. 


richness measures was: 

IV 200 = 3.69±0.0756(AvT’) a669±ool)494 (8) 

Finally we consider the richness and redshift ranges of the 
galaxy clusters, where the redshift used is the photometric redshift 
found using a neural network algorithm. The photo-z used is that 
of the central galaxy in the cluster, the galaxy of highest local den¬ 
sity as determined by the VT cluster finder. For the 19 316 Stripe 
82 co-add clusters the redshift range is 0 < z < 0.98 and the rich¬ 
ness range is 1 < Nyr < 99. For the maxBCG catalog the redshift 
range is 0.1 < z < 0.3 and the richness range is 10 < N 200 < 188. In 
our analysis of the VT clusters, we imposed a criterion that clusters 
needed to be at z = 0.1 or higher (before any photo-z corrections) 
because as redshift approaches 0, the significance of the Nyr de¬ 
tection decreases significantly. In Figures [3] and [4] we present his¬ 
tograms of richness and redshift for the maxBCG and the Stripe 
82 co-add cluster samples. 

Source galaxies were selected from a catalog of galaxies 
found in the Stripe 82 co-add. In order to produce this catalog. 


Weak lensing can cause a systematic change in the ellipticities of 
many galaxies behind the galaxy clusters; this observed change 
is called shear , represented by y. Shear has two components, 
tang ential (y) and orthotangential_or cross shear (y x ), given by 
l iSchneider, Kochanek & Wambgans J2005I) 


y t = —91 ye 2 '0 
Yx = -3 [ye~ 2, ' f 


(9) 


where <j) is the position angle measured between the center of the 
cluster and the lensed galaxy. Only tangential shear is produced 
by weak lensing, so we expect no signal in orthotangential shear. 
Since the observed tangential shear is related to cluster mass, weak 
lensing can provide a measure of cluster mass. Tangential shear is 
found for a combination of galaxy clusters and source galaxies, 
where source galaxies are galaxies nearby clusters in projection, 
but at a larger redshift than the clusters. The source galaxies are 
the galaxies whose shape is measured, and the mass associated 
with the clusters is the cause of the tangential shear that describes 
the shape variation. 

The quantity that can be directly observed in images is galaxy 
ellipticity, which we found by measuring second order moments of 
galaxies. In order to convert this to shear, we divided by respon- 
sivity (R). For the Stripe 82 co-add, we used 

/? = 2(1 — ofjv) = 1.73 (10) 


where o$y is the intrinsi c galaxy sh ape noise in the SDSS, taken 
to be (Jsn = 0.37 jHirata et al.ll2004l) . Shape noise describes the 
intrinsic variation in galaxy shapes, independent of any measure¬ 
ment uncertainty. For Stripe 82 data, PSF deconvolution w as do ne 
by R. Reis using the Hirata-Seljak method (!Hi rata & Seliaki2003h . 

Even in the Stripe 82 co-add data, which is much deeper than 
the SDSS main sample, the weak lensing signal was likely to be 
noisy for an individual cluster. Thus in order to maximize shear 
signal to noise we measured stacked weak lensing shear (following 
Ijohnston et al .1 12007 )). To stack weak lensing shear we measured 
shear for all source galaxies found in a particular richness and red¬ 
shift bin. Thus we took the average of the source galaxy shears in 
that richness/redshift bin at each radius bin. For each bin we had 
a set of 30 measurements of average shear for 30 increments of 
distance from the central BCGs (from 0 — 3.0 h ~ 1 Mpc in steps of 
0 . 1 ). 

Tangential shear is related to the mass and concentration of 
the lensing mass distribution. In order to find these quantities, we 
apply an NFW profile and find the bes t-fitting mass and concen - 
tration that produce the observed shear (IWright & Brainerdl200(j) . 
Weak lensing shear as a function of radius predicted by the NFW 
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c 200 as a function of M 200 



m 200 ( 1 0“ h“' M 0 ) 

Figure 5. A plot of C 200 as a function of M 200 using the relationship de¬ 
scribed in lDuffv et alj 12008I) . The uppeiTnost (black) curve is for median z 
of 0.25, the middle (red) curve is for median z of 0.55 and the lowest (blue) 
curve is for median z of 0.85. 


Tangential shear vs. radius for Low-z, Nvr =2 



radius (Mpc) 

Figure 6. A sample shear profile. This profile is for Stripe 82 co-add, low-z 
bin (0.1 < z < 0.4), Nvt — 2 with no corrections for systematics applied. 
Note that the upper tangential shear profile shows a shear signal while the 
lower orthotangential shear profile does not. There were 614 clusters in 
this bin and 2 699 077 source galaxies. 


model is given as 


Ynfw{x) = 


Znfw (x) — Enfw (•*) 


( 11 ) 


Here x is a dimensionless radius equal to r/r s , where r is the dis¬ 
tance from the cluster center in the lens plane and r s is the scale 
radius in the NFW model, the radius at which the density changes 
from a 1/r relation to a 1 /r 5 relation. T-nfw (x) is the surface mass 
density of the galaxy cluster’s dark matter halo in the NFW model, 
£ NFw(x ) is the mean surface mass density of the halo and £ cn? is 


£<-rif — ■ 


c 2 d s 


( 12 ) 


4 nGdidLs 

Here ds, dL and dis are angular diameter distances to the source 
galaxy, to the lensing cluster and from lensing cluster to source 
galaxy, respectively. We implemented the series of relations in 
Wright & Brainerd (2000) predicting shear as a function of radius 
in the NFW model. We then fit observed shear to this model in or¬ 
der to find the mass QVHoo) thst would lead to the measured values 
of tangential shear in our richness and redshift bins. We held c '200 
constant while doing these fits, as we found that if we allowed 
both M 200 and C 200 to float simultaneously, values for C 200 were 
lower than would be expected and were not well constrained. We 
cho se val ues of £200 by using the concentration-mass relation of 
iDuffv et all d2008l) . given as: 


C200 (M,Z)=A (M/M pivot ) B ( 1 + zf 


(13) 


We use the values given for the full (not only relaxed) cluster sam¬ 
ple and redshift range 0 — 2. For this sample, A = 5.74±0.12, B = 
—0.084±0.006. C = —0.47±0.04 and M pivot = 2x 10 12 lr l M Q . 
We found C200 in the range of 10 12 — 10 15 and then took 

the average value of C 200 across this range of values for Miqq and 
for the median redshift value in each redshift bin (0.25, 0.55, 0.85 
and 0.55 for the low-z, mid-z, high-z and all-z bins, respectively). 

We found C 200 = 3-33 for the low-z bin, C 200 = 3.01 for the 
mid-z and all-z bins and C 200 = 2.77 for the high-z bin. In Figure 
[5]we plot C 200 as a function of M 200 f° r the three median redshifts. 

We then calculated a scale factor to scale all cluster shears to 


the same redshift. This was done because weak lensing shear is a 
function of redshift, and while all clusters in a redshift/richness bin 
are at similar redshift, they are not at identical redshift. The scale 
factor S is 

c _ d sdL.fiddLS_fid (W) 

ds.fiddLdLS 

The ones marked ”fid” mean at the fiducial redshift, which is 
Zfiducial Jens 0.55 for the lens and Zfiducial source 0.75 for the 
source. Thus all objects are scaled to these two redshifts, which 
were chosen as approximately the mean redshifts for the cluster 
and the source galaxies respectively. Equation |T4] was obtained by 
taking the ratio of £ eri , for the cluster redshift to that the fiducial 
redshift. We calculate average tangential shear by using inverse 
variance weighting in order to overcome the effects of large values 
of S that arise when di and d$ are close in value. 

Finally we take the measured average shear values in each 
bin in the range of 0 — 3.0 hr 1 Mpc from the central BCGs and the 
predicted shear values found from Equation QT] and minimize the 
X 2 function relating them. The value of distance from the cluster 
center in each bin was the midpoint of the distance bin; for exam¬ 
ple, in the distance range of 0.1 —0.2 /t -1 Mpc, the distance used 
for fitting was 0.15 h~ x Mpc. The values of M 200 that minimize 
the x 2 are taken as the NFW model fits for cluster mass. A sample 
plot of shear profiles with best-fitting curves is shown in Figure [6] 
for the low redshift (0.1 —0.4) bin with Nyr = 2. Note that there 
is good evidence for a tangential shear signal, but orthotangential 
shear is consistent with zero. 

In Tables |T] and [2] we present measurements of richness and 
stacked cluster mass as well as reduced x? f° r a fit to 0 for each 
richness and redshift bin. The value of Nvt given is the median of 
the Nvt values in that richness bin. The error on Nvt is the range 
of Nvt values, given by 

. . 1 tty t jnax ~ tty I jnin , c 

A Ny T = - - - (15) 

The errors on the M 200 values are standard deviations on values 
output by the fitting routine. The x 2 ec i we report here is the reduced 
%red f° r the nu U test of fitting the weak lensing shear to 0. 
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Table 1. Weak lensing shear fit results for low-z (0.1 —0.4) and mid-z 
(0.4 — 0.7) clusters. This table includes values for richness (Nvt), mass 
(M 200 , in units of 10 14 hr 1 Mg) and null test (fit to 0) reduced % 2 . Values 
for concentration (C 200 ) are 3.33 for low-z and 3.01 for mid-z. No correc¬ 
tions for systematics are applied to obtain these values. 


Nyj 

M 200 

Xred 

LOW-Z 

1 ±0 

0.0169 ±0.0100 

1.57 

2 ± 0 

0.0471 ±0.0135 

1.84 

3 ± 0 

0.0797 ± 0.0200 

2.37 

4 ± 0 

0.0695 ± 0.0238 

1.66 

5 ± 0 

0.131 ±0.0321 

2.32 

6 ± 0 

0.117 ±0.0361 

1.67 

7 ±0 

0.181 ±0.0510 

1.98 

8 ±0 

0.130 ±0.0496 

1.97 

9 ± 0.500 

0.207 ± 0.0478 

2.32 

11 ±0.500 

0.464 ± 0.0773 

4.26 

14 ± 1.00 

0.377 ± 0.0749 

3.91 

18 ±2.00 

0.288 ± 0.0706 

2.42 

23 ± 4.50 

0.421 ± 0.0973 

2.68 

36 ± 4.50 

0.904 ± 0.204 

2.62 

46 ± 4.50 

1.21 ±0.258 

2.49 

54 ± 4.50 

1.92 ±0.532 

2.10 

66 ± 3.00 

1.69 ±0.536 

2.15 

78 ± 4.00 

2.081 ± 0.870 

1.81 

99 ±0 

0.560 ± 0.842 

1.33 

MID-Z 

1 ±0 

0.0614 ± 0.0229 

1.49 

2 ± 0 

0.00693 ± 0.00991 

1.55 

3 ±0 

0.0486 ± 0.0240 

1.23 

4 ± 0 

0.173 ±0.0459 

1.91 

5 ±0 

0.130 ±0.0526 

2.084 

6 ± 0 

0.217 ±0.0708 

2.016 

7 ±0 

0.204 ± 0.0802 

1.44 

8 ±0 

0.506 ±0.132 

1.91 

9 ± 0.500 

0.431 ±0.110 

2.22 

11 ±0.500 

0.654 ±0.166 

2.21 

14 ± 1.00 

0.0498 ± 0.0672 

0.963 

17 ± 2.00 

0.605 ±0.171 

1.93 

24 ± 4.50 

0.799 ± 0.233 

1.57 

33 ± 3.50 

0.922 ± 0.579 

1.57 

44 ± 3.00 

1.042 ±0.817 

1.61 

55 ± 2.00 

1.643 ±0.797 

1.41 

63 ± 4.00 

1.97 ± 1.20 

1.069 

91 ± 7.50 

0.738 ± 1.21 

0.712 


The xt-ej we re P or t here is not the same as the x 2 ,,,i used for 
the fit to the NFW profile. What we are checking here is whether 
the tangential shear profile is well fit by zero shear signal. If this 
is a good fit (reduced x~ ~ 1) then there is no signal. If this is a 
bad fit (reduced x 2 > 1). then we have cause to believe there is a 
tangential shear signal. Note that in the low-z and mid-z bins the 
X 2 ed are mostly well above 1, which suggests that the fit to 0 is 
poor and there is measurable weak lensing signal. For the high-z 
bin many of the x},,,i values are lower, which suggests the weak 
lensing signal is less clear for the high-z bin. 


5 CONSIDERATION OF SYSTEMATIC ERRORS 

We next consider tests to verify our measurements of tangential 
shear. We also consider systematic errors that can bias measure¬ 
ments of shear and thus measurements of mass and concentration 


Table 2. Weak lensing shear fit results for high-z (0.7 — 1.0) and all-z (0.1 — 
1.0) clusters. This table includes values for richness (Nvt), mass (M 200 , 
in units of 10 14 h " 1 Mg) and null test (fit to 0) reduced ■ Values for 
concentration (C 200 ) are 2.77 for high-z and 3.01 for all-z. No corrections 
for systematics are applied to obtain these values. 


Nyj 

M 200 

tCred 

HIGH-Z 

1 ±0 

0.0200 ± 0.0375 

1.01 

2 ± 0 

0.0519 ±0.0505 

0.931 

3 ±0 

0.165 ±0.0949 

1.34 

4 ± 0 

0.0343 ± 0.0612 

1.18 

5 ±0 

0.312 ±0.173 

1.91 

6 ± 0 

0.281 ±0.196 

0.774 

7 ±0 

0.0533 ±0.128 

0.562 

8 ±0 

0.00642 ± 0.0661 

1.01 

9 ± 0.500 

0.680 ± 0.355 

0.983 

11 ±0.500 

0.258 ± 0.325 

0.656 

14 ± 1.00 

0.212 ±0.343 

1.03 

17 ± 2.00 

1.56 ±0.917 

0.886 

24 ± 4.50 

1.66 ± 1.24 

1.34 

34 ± 4.50 

5.31 ±4.49 

1.18 

44 ± 5.00 

0.0483± 1.023 

1.20 

ALL-Z 

1 ±0 

0.0288 ± 0.00986 

1.64 

2 ± 0 

0.0379 ±0.0102 

1.77 

3 ±0 

0.0728 ± 0.0158 

2.96 

4 ± 0 

0.100 ±0.0217 

2.74 

5 ±0 

0.144 ±0.0282 

3.20 

6 ± 0 

0.150 ±0.0330 

2.95 

7 ±0 

0.174 ±0.0417 

2.71 

8 ±0 

0.212 ±0.0502 

2.58 

9 ± 0.500 

0.276 ± 0.0460 

3.46 

11 ±0.500 

0.515 ±0.0716 

5.70 

14 ± 1.00 

0.304 ± 0.0618 

4.052 

17 ± 2.00 

0.390 ± 0.0702 

3.73 

24 ± 4.50 

0.533 ± 0.0950 

3.62 

35 ± 4.50 

0.980 ± 0.202 

2.78 

45 ± 4.50 

1.21 ±0.253 

2.44 

54 ± 4.50 

1.93 ±0.461 

2.60 

63 ± 4.00 

1.82 ±0.511 

2.43 

78 ± 4.00 

2.12 ±0.893 

1.72 

99 ± 4.00 

0.877 ± 0.826 

0.986 


in the NFW model. Our tests to verify our measurements of shear 
include measurement of orthotangential shear and measurement 
of tangential shear around random points rather than nearby clus¬ 
ters. Systematics we considered included photo-z bias, the effect 
of the central BCG, dark matter halo miscentering, uncertainties 
on photo-zs and the contamination of source galaxies by misiden- 
tified lensing galaxies. 

5.1 Orthotangential Shear Test 

For each richness and redshift bin we made plots of both tangen¬ 
tial and orthotangential shear. For tangential shear we expect to 
see an increase in shear at small distances from the cluster center 
(consistent with non-zero tangential shear) while for orthotangen¬ 
tial shear we expect to see no change in the orthotangential shear 
at small distances (consistent with zero orthotangential shear). If 
we see significant signal in the orthotangential shear, this can be 
evidence of error as weak lensing will not produce orthotangential 
shear. As seen in Figure [6] typically we can see a clear signal in 


© 2014 RAS, MNRAS 000.1711141 


Mass Calibration in Stripe 82 7 


Low —z 



0.00 0.02 0.04 0.06 0.08 0.10 0.12 
Tangential 

High-z 



Tangential 


Mid-z 



Tangential 

All —z 



Tangential 


Figure 7. y 2 probabilities compared for tangential and orthotangential 
shear. The green line is the y = x line. 


tangential shear while the orthotangential shear is apparently con¬ 
sistent with zero. 

In Figure [7] we compare the y} probability calculated for the 
tangential shear to that calculated for orthotangential shear. This 
quantity tells us the probability that the y 2 , for a fit to 0 will ex¬ 
ceed the given value of xLj- For null data, like we expect orthotan¬ 
gential shear signal to be, this probability should be on the order 
of 0.5, since we expect the orthotangential shear to be well fit to 0. 
For tangential shear this probability should be very small, since we 
would expect tangential shear measurements to not be consistent 
with 0. Note that in Figure[7]we find that y 2 probability is almost 
always significantly higher for orthotangential shear than it is for 
tangential shear. For orthotangential shear, the probability ranges 
from 0 to 1, and so on average is about 0.5, as expected. In most 
cases, the probabilities for tangential shear are near 0, meaning 
they are not well fit to a tangential shear of 0. Thus we conclude 
that in most cases, there is evidence for tangential shear signal, and 
not for orthotangential shear signal. The one notable exception is 
the high-z (0.7 < z < 1.0) data. The high-z data has the lowest 
statistics of any of the bins, and so the tangential shear fits have 
higher error bars. The tangential shear signal is less obvious, with 
y 2 probabilities from 0 to 1, while the orthotangential shear signal 
is unchanged. 


5.2 Random Points Test 

We would expect no tangential shear signal around random points 
in the sky while we would expect a non-zero tangential shear 
around clusters. If random points were to give signal similar to 
that measured near clusters, this would be evidence that there is a 
systematic problem in our measurements. In order to conduct this 
random points test we first generated a sample of random points. 
We generated 19 706 random values for RA and DEC within the 
ranges of the real values, with a distribution of Nyr and z similar 
to that found in the real data. 

Once we generated the random points catalog, we remeasured 
tangential shear using the random points as the clusters. We found 
that the tangential shear profiles generated were now consistent 
with zero, as would be expected for random points. In Figure [8] 
we plot the reduced y 2 values from tangential shear fits to zero 
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Figure 8. Reduced yy values for tangential shear fits to zero for clusters 
and for the same fits to random points, given compared to Nvt- Note that 
for the clusters, the reduced yy values are well above 1, meaning they are 
not consistent with a model of no shear signal. For the random points, 
reduced yy values are close to 1, suggesting no evidence of shear signal 
(as expected). 


for both the actual data (upper) and for the random points (lower). 
Note that for the real data the reduced y 2 values are mostly larger 
than 1, meaning that the real data is not consistent with fit to zero, 
that is a model of no shear signal. On the other hand, for the ran¬ 
dom points the reduced y 2 values are mostly near 1. meaning that 
the random points are consistent with no tangential shear. This re¬ 
inforces our claim that we have measured a significant weak lens- 
ing shear signal. 


5.3 Photo-z Bias 

A photometric redshift (or photo-z) is a redshift measured by ob¬ 
serving magnitudes in multiple filters and then using a neural net¬ 
work or other method to estimate redshift based on this informa¬ 
tion and a training set of spectroscopic redshifts (spec-z). The re¬ 
lation between the magnitudes measured in different filters has a 
dependence on the redshift of an object, thus photometry can be 
used to estimate redshifts. Photometric redshifts are usually used 
in large surveys where it is implausible to find spectroscopic red- 
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Figure 10. M 200 with and without the central 0.1 h 1 Mpc included in 
fitting. The green line is the y = x line. Note that the points mostly follow 
this line. 


Figure 9. A comparison of photometric redshift to spectroscopic redshift 
for clusters that had at least one member with a spectroscopic redshift. The 
photo-zs were the redshift of the central BCG in each cluster. This relation 
was used to find photo-z bias corrections for each redshift bin. 


shifts for most objects. The advantage to this method is that red- 
shifts can be estimated for many objects, but the disadvantage is 
that errors on photometric redshifts are significantly larger than 
those for spectroscopic redshifts. We consider two systematic ef¬ 
fects caused by photo-zs. In this section we consider a small bias 
in the photo-zs and in Section l5^6l we consider the effect of photo-z 
uncertainties on the value of £~1 f . 

When comparing the photometric redshifts of VT clusters to 
maxBCG cluster redshifts and spectroscopic redshifts for some 
cluster galaxies we found a small bias (see Figure [5). We found 
that the cluster photometric redshifts were too high by 0.0277 in 
the low-z bin, by 0.0348 in the mid-z bin, by 0.00974 in the high- 
z bin and by 0.0308 in the all-z bin. In order to correct for this, 
we subtracted these values from all cluster redshifts before these 
redshifts were used for binning and for calculation of E7.L. 

5.4 Central BCG 

In our fitting of the mass and concentration, we fit an NFW pro¬ 
file to the dark matter halo of the clusters. In so doing, we did not 
include a term for the mass of the central BCG. Thus we wanted 
to try to estimate what effect the mass of the BCG would have on 
the shear profile. To do this we excluded the central 0.1 Mpc. 
We did this as the mass of the BCG would usually affect only the 
innermost regions of the shear profile. In Figure [TO] we compare 
the cluster masses found with and without the central 0.1 /j _1 Mpc 
region of the clusters. The green line is the y = x line. The points 
follow this line quite closely, with a few exceptions. Thus we con¬ 
clude that the mass of the central BCG has minimal impact on 
the shear profiles. Nevertheless for our final measurements we ex¬ 
clude the central 0.1 Mpc of the stacked cluster profiles from 
the NFW fits in order to avoid bias from the mass of the central 
BCG. 


5.5 Halo Miscentering 

When we fit an NFW profile to the dark matter halos of the galaxy 
clusters, we take the BCG to be the center of the halo, measur¬ 
ing all distances with respect to this. However, as mentioned in 
Section [2] this will be incorrect some fraction of the time, as the 
BCG at times is not the center of the dark matter halo. If mass is 
measured at a location outside the center of a halo, the mass found 
for that halo will be less than it should be, since the outer region 
of a halo is then being taken as the center. As this will bias mass 
measurements, our analysis needs to include a consideration of the 
contribution from halos that are not centered on the BCG. In or¬ 
der to consider halo miscentering we use Equation [3l lSimet et al.l 

I201 A _ 

ISimet et al. J2012h base their method on the maxBCG 
method, and so their probability distribution of miscentering 
in dark matter halos will follow the distributions discussed in 
Ijohnston et al. < 2007 ). We examined whether the probability of 
miscentering of clusters in halos would be the same for our VT 
method as it was for the maxBCG method. We found a very simi¬ 
lar miscentering_distributionusing simulations described and used 
in ISoa res-Santo s et ahl 1201 ill , with a 2D Gaussian fit giving a 
standard deviation a = 0.47 /i _1 Mpc while ljohnston et all J2007h 
found a = 0.42 /j~* Mpc. 

We applied Equation [3] by taking the masses output by the 
fitting routine as M 200 ,mts and multiplying them by 1.44. Then the 
true masses (Miooyrue) are a function of A^qo, following Equation 
[3] Before we calculated M2oo,true we had to convert Nyr to A 200 
using Equation [8] In Figure QT| we compare masses found without 
the miscentering correction to masses found after the miscentering 
correction. Note that all masses are systematically higher. For each 
of the redshift bins we found the ratio of miscentering-corrected 
mass to the uncorrected masses and then found the median of 
the ratios. This median is shown in Figure |TT] as the quantity in¬ 
crease. Note that the percent increases are signficant and range 
from 46.6% to 53.3%. 
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Figure 11. M 200 values before and after miscentering correction. The green 
line is the y — x line. Note that all masses are pushed upward by about 50%. 


5.6 Uncertainties on Photometric Redshifts 

All source galaxies in the sample had photometric redshifts mea¬ 
sured for them. Our weak lensing mass measurements depend on 
photometric redshifts because tangential shear / depends on If/,, 
which is a function of redshift. Thus we need to consider the errors 
which may be introduced into our measurements by uncertainties 
on photometric redshifts. To do this, we assembled a sample of 
spectroscopic redshifts to compare to. 

We gathered public data from four data sets: The Stripe 82 
co-add spectroscopic data, VVDS, DEEP2 and VIPERS. To ac¬ 
cess spectroscopy from the Stripe 82 co-add we used the SDSS-III 
catalog archive server to query for spectroscopic redshifts for ob¬ 
jects within Stripe 82 that had i-band magnitudes 19 < i < 23.5, 
type=3 (galaxies) and zwarning=0. The last paramete r is a mea sure 
of the quality of the spectroscopic redshift. DEEP2 dNewman et all 
120131) is a project which uses the Keck telescopes at Mauna Kea. 
The DEEP2 project seeks to understand the evolution of galaxies 
and galaxy clusters by taking spectra of more than 50 000 galax¬ 
ies out to red shift grea ter than 0.7. The VIMOS VLT Deep Sur¬ 
vey (VVDS) dLe Fevre et al.M2013l ) is a survey conducted using 
the VIMOS (Visible imaging Multi-Object Spectrograph) instru¬ 
ment on the Very Large Telescope (VLT). VVDS observed 3 fields 
that observed 35 526 galaxies of magnitude ranging between 17.5 
and 24.75 and redshift between 0 and 6.7. The VIMOS Public 
Extragalactic Redshift Survey (VIPERS) dOarilli et~aPl2014l) is a 
survey of high-redshift objects taken using the VIMOS spectro¬ 
graph on the VLT. VIPERS observed two fields of objects cho¬ 
sen from the Canada-France-Hawaii Telescope Legacy Survey- 
Wide (CFHTLS-Wide) optical photometric catalog. In the first 
public data release of VIPERS there are 54 756 galaxies with 
spectroscopy. In each of these surveys we match the galaxy cat¬ 
alog with our own catalog of galaxies in Stripe 82. Our final cat¬ 
alog of objects with spectroscopic redshifts is a list of all objects 
from these 4 sets of data that overlap with objects in the Stripe 82 
co-add. We found a total of 34 033 matched source galaxies with 
spectroscopic redshifts. 

In order to better constrain the effect of photo-z error, we an¬ 
alyzed how E“], changed with different measures of redshift. We 
looked at four different ways of finding Eli and then we com¬ 
pared them to see how much E/, varied based on the method: 


(i) Find E i using the median photo-z of the source galaxies. 

(ii) Find E“J, for each source galaxy and then take the average, using 
photo-zs. 

(iii) Find E“i for each galaxy that has a spec-z but use the photo-z 
for each of those galaxies. Photo-zs must be weighted to match 
the distribution of the much larger Stripe 82 photo-z sample (see 
details below). 

(iv) Find Eli for each galaxy that has a spec-z but use the spec-z for 
each of those galaxies. Spec-zs must also be weighted, as in (iii). 


We would expect method (iv) to give the best results since it 
uses spectroscopic redshifts, but since the number of spectroscopic 
redshifts is small, we have to first demonstrate that the weighting 
method works. In order to demonstrate this, we use methods (ii) 
and (iii), which both use photo-zs. Method (ii) involves finding 
E/ for each source galaxy and method (iii) involves finding it for 
galaxies that have spectroscopic redshifts but with weighting. If 
they match, then this suggests that the weighting works well, and 
can then be used with the spectroscopic redshift sample. 

The original method (method (i)) by which we found E~j, 
was to find the photometric redshift of the cluster and then find 
the median of all of the photometric redshifts of the source galax¬ 
ies. These redshifts were used to find the angular diameter dis¬ 
tances, then E7.1 was calculated from them. Galaxies were con¬ 
sidered source galaxies if they were (1) within 3 h~ l Mpc of the 
BCG, (2) had i-band magnitude between 20.5 and 23.5, (3) had 
ellipticities less than 2 and (4) had Zsource > Zduster + 0.1. To test 
this method we produced a table of 1000 redshifts from 0.1 — 1.0, 
treating these as cluster redshifts. We then took the first 1000 ob¬ 
jects in the source galaxy catalog and looped through these to find 
source galaxies to correspond to each of the fake clusters; we only 
required that the source galaxies have 20.5 < i < 23.5 and that 
Zsource > Zduster + 0.1 • We then took the median photo-z of these 
source galaxies and used this number along with redshift of the 
cluster to find E~ L. 

For method (ii) (using all photometric redshifts) we used al¬ 
most the same method as in the previous Section. However rather 
than finding the median source redshift, for each cluster we looped 
through all source galaxies and found EtI for each. Then for each 
cluster we found a mean value of E“L. 

In method (iii) we began to use the 34 033 galaxies for which 
we had both photometric and spectroscopic redshifts. However 
since we were using this small spectroscopic sample to find val¬ 
ues of E“), for the full Stripe 82 co-add dataset, we had to weight 
the small spectroscopic sample to match the distribution of the 
much larger Stripe 82 sample. We utilized a routine in C++ called 
CALCWEIGHTS.CPP developed by C. Cunha that uses a nearest- 
neighbor code to extrapolate a redshift distribution ( iLima et all 
120081 ) . 

Now that we had these weights we could find E“j, using 
method (iii). We produced a cluster sample like before, ranging 
in redshift from 0.1 — 1.0. The source redshifts we used in method 
(iii) were the photometric redshifts that corresponded to the 34 
033 galaxies with spectroscopic redshifts. Then for each cluster 
redshift we looped through all the (photometric) source redshifts, 
producing a value of E er ; f x weight for each. Finally we found the 
mean by finding for N source galaxies: 


/ N \ 

L ^critj x weighti 


N 

E weighti 

i—t 


(16) 
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£ as a function of redshift 



Redshift 

Figure 12. XT" !, as a function of redshift for each of the four different meth¬ 
ods described in this section. The dash-dot violet line is from method (i), 
the solid black line is from method (ii), the dotted red line is from method 
(iii) and the dashed blue line is from method (iv). All values of Z~| ( are 
multiplied by 1CT 16 m 2 /kg. 


Finally in method (iv), we used the same method as method 
(iii), but this time the source redshifts were the spectroscopic red- 
shifts. Again we found x weight for each source galaxy for 
a particular cluster redshift and again we found the mean value of 
Z7.L for that cluster redshift by using Equationll6l 

Our final product in each case was a lookup table of values of 
Z ~ jt for a set of cluster redshifts ranging from 0.1 — 1.0. In Fig¬ 
ure [12] we plot Z“L as a function of redshift from 0.1 — 1.0. In 
this plot the dash-dot violet line is from method (i), the solid black 
line is from method (ii), the dotted red line is from method (iii) 
and the dashed blue line is from method (iv). Note that the origi¬ 
nal method of finding Z~.J f (method (i)) deviates the most from the 
other methods. The red and the black lines ((methods (ii) and (iii)) 
match closely, indicating that the weighting method effectively re¬ 
produces the larger sample. 

We finally conclude that there are significant differences in 
Zdepending on which method of finding it we use. We also 
conclude, from the closeness of the black and the red lines, that 
the weights code is highly effective at reproducing the photomet¬ 
ric distribution from the spectroscopic distribution. Therefore we 
finally chose to use the weighted distribution of spectroscopic red¬ 
shifts (the blue line from method (iv) in Figure 1121 to calculate 
Z~’ 

We produced a lookup table of values of Z~.J f using these 
spectroscopic redshifts and then had the fitting routine refer to this 
table when finding Z“l ( for a particular cluster. The scale factor 
fEquation llTI l was then calculated as: 

^ ^-'crit .actual (17) 

7 crit-fiducial 

Z c „y .fiducial was taken as 7.94 x 10 15 m 2 kg~ l , which is the value 
of L crit found for z SO urce = 0.75 and z c iuster = 0.55. 

5.7 Foreground Galaxy Contamination 

Weak lensing shear was measured for a sample of galaxies that 
were more distant than each of the lensing clusters. We measured 
shear only for galaxies that were at a redshift that was 0.1 greater 


than cluster redshift. The redshifts used for these cuts are pho¬ 
tometric redshifts, and thus they have larger errors than if they 
were spectroscopic redshifts. These errors can allow some galax¬ 
ies that are actually at lower redshift than the cluster (i.e., in front 
of the cluster) to be included in the sample of background source 
galaxies. This would contaminate the shear signal with galaxies 
that cannot exhibit shear. This is especially a problem nearby a 
galaxy cluster, as cluster galaxies can be misidentified as back¬ 
ground galaxies. The contamination can be quantified by noting 
that average tangential shear (y m , ? ) in a particular richness bin is 


A ’real Nfake 

L YTjreal T L Yl-fake 

^41 (18) 

where real and fake refer to source galaxies and misidentified fore¬ 
ground galaxies respectively and A means total number of galaxies 
found in that richness bin. We expect ’LYr.fake t° sum to zero since 
this includes only shape noise and no actual shear. However, the 
Nfake term remains, and we need to find a correction factor C(r) 
to account for it. 


What we want is to remove the Nf a k e term in the denominator. 
Thus we need to multiply the measured y ave by this term: 


C(r) 


threat T ^fake 
threat 


(19) 


But since we cannot directly measure N rea i, we can instead mea¬ 
sure a number density, n, the number of galaxies per area per clus¬ 
ter. Thus 


C(r) 


Nreai T Nfake 
Nreai 


ktmeasuredy) 


^random 


( 20 ) 


Here n measure d(r) is the number of galaxies measured per area 
per cluster at various radii around clusters; this tells us about 
Nreai + Nf a k e because it accounts for all apparently source galax¬ 
ies found nearby clusters (including misidentified cluster galax¬ 
ies). The second term, n ranc i om , is the number of galaxies per area 
per cluster measured at random points; that is, far enough away 
from the center of the cluster that the concentration of galaxies 
that marks a galaxy cluster does not affect the number density. 

Number density n is measured by finding the total number of 
galaxies per area per cluster in a given richness bin. To find this, we 
find the total number of source galaxies in bins of 0.1 /A 1 Mpc of 
distance from the BCG in each richness and redshift bin. We then 
divide this number by the product of the area of an annulus at that 
increment of distance from the BCG and the number of clusters in 
that richness bin. In other words, we measure 


Aqnnu/us 

Nclusters 


( 21 ) 


where N tota i is the total number of apparently source galaxies in 
that richness and redshift bin, A annu j us is the area of the annulus 
and N c i uster s is the number of clusters in that bin. 

When we first measured the number densities, we found that 
they all decreased with radius, even for the random points. This 
is unexpected behavior, as random points should overall have no 
dependence of number density on position. It was found that this 
behavior was due to the nature of Stripe 82 as a thin strip of ob¬ 
servations (2.5 degrees wide in declination). For clusters that were 
near the edge of the images the number of nearby galaxies that 
could be measured would be limited to those that appeared on the 
images. As the distance from the cluster center increased, the num- 
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Figure 13. Plots showing number density as a function of radius found 
for the first six richness bins for 0.4 < z < 0.7. The horizontal green line 
represents the median value of number density and the vertical red line 
represents the beginning radius at which the random number density was 
measured. 


ber of galaxies per area would drop, not due to an actual decrease 
but simply because we had progressed outside the available imag¬ 
ing region. To address this issue, we remeasured the number of 
source galaxies per area per cluster while rejecting any clusters for 
which a circle of radius 3.0 /t~' Mpc would be off the imaging 
area. 

We found that number density for random points in the sky 
is about the same as the number density in the region of a cluster 
but far away from the center. Because of this we chose to use the 
number density in the region of clusters but far away as the refer¬ 
ence point for our foreground corrections. As the number density 
approached a constant by 1.4 h ~ 1 Mpc, we measured the random 
number density between 1.4 — 3.0 Mpc (see Figure fl3l>. 

We finally multiplied the average tangential shear (y fl ,, e ) 
by the C(r) at each radius. The values of C(r) usually ranged 
from about 1.5 at the most to 0.9 at the least. The corrections are 
applied for each radius bin (of size 0.1 /t _1 Mpc), but beyond the 
first several radius bins most corrections are close to 1 (meaning 
no change). 


Table 3. Mass-richness relation fit parameters for the Stripe 82 co-add. 
Mass coefficients all have the units 10 13 h ! Mq. The Xt?d describes the 
goodness of the power law fit to the mass-richness data. 


Redshift Bin 

Mass Coefficient 

Slope of Power 
Law 

Xli 

No corrections 

0.1-0.4 

5.23 ± 0.375 

1.05 ±0.0632 

1.10 

0.4-0.7 

7.46 ±0.917 

0.910 ± 0.0996 

1.28 

0.7-1.0 

16.9 ± 5.75 

1.38 ±0.271 

0.351 

0.1-1.0 

5.78 ±0.368 

1.04 ±0.0548 

1.22 

Photo-z bias 

0 . 1 - 0.4 

4.71 ±0.322 

1.02 ±0.0598 

1.52 

0 . 4 - 0.7 

7.48 ± 1.06 

1.10 ± 0.118 

0.728 

0 . 7 - 1.0 

20.0 ± 6.63 

1.45 ±0.280 

0.548 

0.1-1.0 

5.10 ±0.322 

1.04 ±0.0528 

1.45 

Central BCG 

0 . 1 - 0.4 

5.85 ± 0.402 

1.03 ±0.0642 

1.18 

0 . 4 - 0.7 

8.17 ±0.969 

0.899 ±0.0931 

1.32 

0 . 7 - 1.0 

18.8 ±6.10 

1.31 ±0.259 

0.425 

0.1-1.0 

6.39 ± 0.394 

1.02 ±0.0539 

1.31 

Halo Miscentering 

0 . 1 - 0.4 

7.02 ± 0.495 

0.904 ± 0.0628 

1.13 

0.4-0.7 

10.0 ± 1.22 

0.768 ± 0.0990 

1.28 

0.7-1.0 

22.8 ± 7.73 

1.25 ±0.271 

0.352 

0.1-1.0 

7.76 ± 0.487 

0.901 ± 0.0538 

1.26 

Photo-z uncertainty 

0 . 1 - 0.4 

5.73 ±0.410 

1.06 ±0.0641 

1.10 

0.4-0.7 

8.14 ±0.999 

0.933 ± 0.0998 

1.27 

0.7-1.0 

12.9 ± 4.39 

1.43 ±0.281 

0.347 

0.1-1.0 

6.26 ± 0.401 

1.07 ±0.0547 

1.20 

Foreground galaxies 

0 . 1 - 0.4 

5.66 ±0.391 

1.08 ±0.0612 

1.15 

0.4-0.7 

8.22 ± 0.961 

0.897 ± 0.0948 

1.20 

0.7-1.0 

19.6 ±6.10 

1.33 ±0.240 

0.529 

0.1-1.0 

6.28 ± 0.388 

1.04 ±0.0511 

1.33 

All corrections 

0 . 1 - 0.4 

8.49 ± 0.526 

0.905 ± 0.0585 

1.89 

0 . 4 - 0.7 

14.1 ± 1.78 

0.948 ±0.100 

0.868 

0 . 7 - 1.0 

30.2 ± 8.74 

1.33 ±0.260 

0.674 

0.1-1.0 

9.23 ± 0.525 

0.883 ± 0.0500 

1.94 


6 MASS-RICHNESS RELATIONS 
6.1 Measuring a Mass-Richness Relation 


In order to find a mass-richness relation by using the weak lensing 
shear results, we fit a power relation to M 2 oo and Nvt values. In 
each case we found mass-richness relations for many richness bins 
in four separate cluster redshift bins, 0.1 — 0.4 (low-z), 0.4 — 0.7 
(mid-z), 0.7 — 1.0 (high-z) and 0.1 — 1.0 (all-z). Note that as men¬ 
tioned in Section [2] all mass-richness relations are given in the 
form: 


^200 — M 200 120 


{%) 


( 22 ) 


where a is the slope of the power law. This parameter tells us 
how quickly mass increases with an increase in richness. M 20 o| 2 o 
is the mass coefficient; we define this to be the value of M 2 00 when 
Nvt = 20 . 


Our initial mass-richness relations for the Stripe 82 co-add 
were found with no corrections for the systematics discussed in 
Section [5] We then also found the mass-richness relations for the 
Stripe 82 co-add including corrections for all of the above system¬ 
atics. We corrected for the bias in the photometric redshifts. We ex¬ 
cluded the central 0.1 /z 1 Mpc in the fits, using only source galax¬ 
ies between 0.1 — 3.0 h~ x Mpc from the BCG. We included the 
correction for halo miscentering. We also used the spectroscopic 
redshift lookup table to find values for £ cri/ and we applied the 
foreground galaxy correction factors when calculating y ave . Our 
final plot of mass vs. richness is shown in Figure [14] and the final 
numerical values for the mass-richness relations are shown in Ta¬ 
blets] We note that the errors are smallest for the low-z and all-z 
bins, while they are largest for the high-z bin. 

Finally we present all mass-richness relations found includ¬ 
ing all systematics for (1) low-z, (2) mid-z, (3) high-z and (4) all-z: 
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Low z: M 200 =8.49±0.526 xIO’ 3 (N w /20)°' 905±0 '° 58 ^id z: M 20O = 1 4.' ± 1.78 xIO” (N w /20)° 948 
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Figure 14. The mass-richness relation (M 200 vs. Nyr) for low-z (upper 
left), mid-z (upper right), high z (lower left) and all-z (lower right), in¬ 
cluding corrections for photo-z bias, the effect of the central BCG, halo 
miscentering, photo-z uncertainty and foreground galaxy contamination. 
The dotted red line is the fit for all redshift bins (all-z) and is included for 
comparison in each case. 



1 10 100 

Richness (N„) 





, /„ \ 0.905±0.0585 

(1) M 200 = (8.49±0.526) x 10 13 r'M© 

,, , /„ /o.948±0.100 

(2) M 200 = (14.1 ±1.78) x 10 13 /r'Mgte 

,, , (l.33±0.260 

(3) M 200 = (30.2±8.74) x 10 13 h~ l M Q 

,, . X 0.883±0.0500 

(4) M 2Q0 = (9.23 ±0.525) x 10 13 h~ l M Q J 

(23) 
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6.2 Comparison of Results 

6.2.1 The Effects of the Systematics 

In Figure |T5] we compare all of the values of the mass coefficient 
and the power law slope from Table [3] The situations in these 
plots are: (1) no corrections; (2) photo-z bias correction; (3) central 
BCG correction; (4) halo miscentering correction; (5) photo-z un¬ 
certainty correction; (6) foreground galaxy correction; and (7) all 
corrections. The correction that typically has the most significant 
effect on the mass coefficient and the power law slope is the cor¬ 
rection for halo miscentering. Mass coefficient also changes sig¬ 
nificantly when the five corrections are all applied together. 


6.2.2 Comparison to Other Groups 

Directly comparing mass coefficients and power law slopes to 
other groups is difficult because other samples have a number of 
elements that mak e them differen t . Our meas urements are differ¬ 
ent from those of Ijohnston et al,l 120071) and Isimet et all d2012t) 
because they use the richness measure N 2 qq while we use Nyr- 
The Voronoi tessellation cluster sample will have different proper¬ 
ties than the maxBCG sample used by these groups as the former 
is independent of galaxy color while the latter relies on it. Sim¬ 
ilarly th e 3D-Matched-Filte r (3 D-MF) clus t er-findi ng algorithm 
used by iFord et all d2014al ) and [Ford et all d2014bl) is indepen¬ 
dent of galaxy color, but its method of searching for galaxy clus- 


Figure 15. Comparison of Stripe 82 mass-richness parameters for all dif¬ 
ferent systematic corrections. The upper four plots compare the mass co¬ 
efficients and the lower four plots compare the power law slopes. The sit¬ 
uations are: l=No corrections; 2=Photo-z bias correction; 3=BCG correc¬ 
tion; 4=Halo miscentering correction; 5=Photo-z uncertainty correction; 
6 =Foreground galaxy correction; 7=A11 corrections. Note that the final (all 
corrections) result from this paper is the large blue circle shown as situation 
7. 


ters based on a fiducial luminosity profile will produce a sample 
with different properties than the VT cluster finder. For exam¬ 
ple, the 3D-MF cluster finder produces a cluster sample with a 
much larger richness range than the VT cluster finder. IFord et all 
d2014bl) discuss the fact that their mass coefficient is quite low 
compared to other measurements; they consider this to be be¬ 
cause their cluster finder locates more low-mass clusters than other 
cluster finders or because it d etects m ore faint galaxies, leading 
to hig her richness val u es. iReves et aid d2008l) use the same sam¬ 
ple as Ijohnston et aid d2007 ) but they define M 2 oo based on p 
(mean density of the Univer s e) rathe r th an p nr j, (critical d ensity 
of the Universe). ( Ford et aD l2014b ) and Isimet et aid d2012l) use 
critical density in their definitions of M 2 qq.). Since mean den¬ 
sity is slightly larger than critical density, we expect that the 
mass values found by Reyes would be slightly larger; IReves et al.l 
(2008) predict that m ass values found usin g mean density will 
be 30% higher. Since Ijohnston et ah| 12007) give a mass coeffi¬ 
cient of 0.88 x 10 14 h Ll M e and IReves et al.l J2008ll give one of 
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1.42 x 10 14 Mq (61% difference), these values approximately fol¬ 
low the expected trend. 

In order to make an approximate comparison, we use Equa¬ 
tion [8] to find that when N 200 = 20, Nyr = 12.5. Thus we find 
mass coefficients at Nyj = 12.5, approximating that these cor¬ 
respond to (V 200 = 20. For the low-z bin we find ^ 2001 ^ 200=20 = 
5.35 ±0.299 x 10 13 This is smaller than the previous 

measurements, including 2hf 2 oo|20 = 8.8± 1.2 x 10 13 /r'Mg from 
Johnston et al., Af 2 oo |20 = 9.56±0.75 x 10 13 from Simet 

et al andM 2 oo |20 = 14.2±0.88 x 10 13 Mq from Reyes et al. For the 
mid-z bin we found M 2 oo|Ar 200 = 20 = 8.6 8 ±0 .83 5 x 1 0 13 
In this approximate redshift range iFord et all d2014al fbh and found 
the significantly smaller values of M 200120 = 2.2 ±0.2 x 10 13 Mq 
or ATinnpn = 2.7+Q4 x 10 13 Mq. However, as mentioned above, 
IFord et all ( 12014al ibi) note that their values of mass coefficient are 
lower than others. 

If we directly compare our values of power law slope to those 
of other groups, we find that our value for the redshift 0.1 — 0 .4 
bin of 0.90 5 ±0.0585 i s close to the values found by Simet et al. 
d2012l) and iReves et all d2008l) of 1.10±0 . 12 an d 1.16±0.09. It 
is farther from the value of Ijohnston et al.1 d2007t) of 1.28 ±0.04. 
Our values of power law slope for the redshift 0.1 — 1.0 bin of 
0.883 ± 0.0500 are lower than those found bv lFord et al.l d2014al) 
and lFord et alJd2014bl) . which were 1.5 ±0.1 and 1.4 ±0.1. 


7 REDSHIFT EVOLUTION 

To examine the evolution of the mass-richness relation, we consid¬ 
ered the mass coefficient and the power law slope as a function of 
redshift. In Figure |T6] we plot the mass coefficient and the power 
law slope vs. redshift for three redshift bins: low-z (2 = 0.25), 
mid-z (2 = 0.55) and high-z (2 = 0.85). We then found linear rela¬ 
tions between mass coefficient (M2qo\n vt= 2o) an d z an d power law 
slope, (a) and z as: 

M 200 \n vt =20(z) = (21.3± 1.66)z±(3.13±5.74) 

(24) 

a(z) = (0.382±0.117)z± (0.799±0.307) 

Figure |T6] suggests that there is strong evidence of redshift evolu¬ 
tion in the mass coefficient. The high-z point has large error bars, 
but the low-z and mid-z points have smaller error bars and also 
give evidence of redshift evolution, leading the overall error on 
the slope of the mass coefficient-redshift relation (Equation 124b 
to be less than 10% of the slope . Our finding here disagrees with 
the finding of IFord et alj d2014bl) : Figure 8 in IFord etafl d2014bt) 
shows the opposite trend, of an overall decrease in mass coefficient 
with redshift. However their error bars are large and they state that 
their negative slope is consistent with zero. We find here that red¬ 
shift evolution of the power law slope is less clear. While there 
is some evidence of redshift evolution of the power law slope, the 
error on the slope of the second relation in Equation [24] is about 
30%. 

In Figure |T6] it does appear that there is redshift evolution 
of the mass coefficient and there may be redshift evolution of 
the power law slope. This would suggest that clusters of the 
same richness are more massive at higher redshift and that at 
higher redshift, cluster mass increases faster with increasing 
richness. This would be consistent with the idea that dark matter 
haloes associated with galaxy clusters contain more dark matter 
at higher redshift. However it is important to note that there 
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Figure 16. The mass at Nvt — 20 as a function of redshift (upper plot) 
and the power law slope as a function of redshift (lower plot) for the Stripe 
82 co-add. The green lines are the linear best-fitting relations, with their 
respective equations displayed at the top of each plot. 

is a difference between real mass evolution with redshift and 
apparent mass evolution caused by measurement systematics. 
One particular systematic that might cause the appearance of 
redshift evolution is the limiting magnitude for cluster finding. 
The Voronoi tessellation method uses an apparent magnitude cut 
rather than an absolute magnitude cut. This means that the sample 
of clusters will be less complete at higher redshift, where galaxies 
will be biased toward dimmer apparent magnitude simply because 
of their large distances. In our sample there are more clusters 
found at higher redshift than at lower redshift because there is a 
larger volume sampled at higher redshift. However at increasing 
redshift, we expect the detection completeness of the cluster 
sample to decrease. Thus there would be poorer statistics at higher 
redshift, which could bias the ultimate measurement of mass 
coefficient at higher-z. In order to further test this mass-redshift 
relation, it would be better to also find this relation using a cluster 
sample that is limited by absolute magnitude rather than apparent 
magnitude. 
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8 CONCLUSIONS 

We present mass-richness relations for the Sloan Digital Sky Sur¬ 
vey Stripe 82 co-add. These mass-richness relations are presented 
for four redshift bins, 0.1 < z < 0.4, 0.4 < z < 0.7, 0.7 < z < 1.0 
and 0.1 < z < 1.0. We present a sample of galaxy clusters and 
we describe how these clusters were found using a Voronoi tessel¬ 
lation cluster finder; cluster richness is given in terms of a rich¬ 
ness measure we call Nyj. Our sample includes 19,316 clusters 
in Stripe 82 between 0.1 < z < 0.98 and 1 < Nyr < 99. We de¬ 
scribe how we measured stacked weak lensing shear in a sequence 
of redshift and richness bins. We found values for M 200 by fitting 
measured shear to an NFW model. 

We describe tests conducted to verify the strength of the weak 
lensing signal; these included measurements of orthotangential 
shear for all cluster bins and measurements of weak lensing shear 
around random points. We also describe the effects of systematics 
on weak lensing shear and thus cluster mass (M 200 ) results. These 
systematics include a bias in photometric redshifts, the effect of the 
central BCG, halo miscentering, photometric redshift uncertainty 
and foreground galaxy contamination. We present methods to cor¬ 
rect for each of these systematics while measuring A^oo- We find 
that the halo miscentering correction has the largest effect on mass 
measurements while the other corrections have smaller effects. 

We present values for the mass coefficient and the power 
law slope for mass-richness relations found using power law fits 
to Stripe 82 data for each of the four redshift bins. After con¬ 
sidering all corrections we find values of the mass coefficient 
(M 20 o|2o) °f 8.49 ±0.526, 14.1 ± 1.78, 30.2 ±8.74 and 9.23 ± 
0.525 x 10 13 /j-'Mq for each of the four redshift bins respectively. 
We find values of the power law slope (a) of 0.905 ± 0.0585, 
0.948±0.100, 1.33±0.260 and 0.883±0.0500 respectively. We 
also considered mass evolution of the mass-richness relation. We 
conclude that there is good evidence of redshift evolution of the 
mass coefficient and less clear evidence of redshift evolution of the 
power law slope. However this apparent redshift evolution may be 
a cluster selection effect. 
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